Equation of state of cubic boron nitride at high pressures and temperatures 
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We report accurate measurements of the equation of state (EOS) of cubic boron nitride by x-ray 
diffraction up to f 60 GPa at 295 K and 80 GPa in the range 500-900 K. Experiments were performed 
on single-crystals embedded in a quasi-hydrostatic pressure medium (helium or neon). Comparison 
between the present EOS data at 295 K and literature allows us to critically review the recent 
calibrations of the ruby standard. The full P-V-T data set can be represented by a Mie-Griineisen 
model, which enables us to extract all relevant thermodynamic parameters: bulk modulus and its 
first pressure-derivative, thermal expansion coefficient, thermal Griineisen parameter and its volume 
dependence. This equation of state is used to determine the isothermal Griineisen mode parameter 
of the Raman TO band. A new formulation of the pressure scale based on this Raman mode, using 
physically-constrained parameters, is deduced. 



PACS numbers: 62.50.+p,61.10.Nz,65.40.De 

I. INTRODUCTION 

Cubic boron nitride (c-BN) is a material with remark- 
able properties: extreme hardness, chemical inertness, 
large band gap and high mechanical and thermal stabili- 
ties. This makes it very attractive for a number of appli- 
cations including abrasive or protective-coating material 
and microelectronic devices. Knowing the properties of 
c-BN under conditions of high pressure (P) and temper- 
ature (T) is important for some of these applications. It 
has also been recognizedA^ that c-BN could be a useful 
pressure sensor in high pressure (HP)- high temperature 
(HT) experiments in diamond anvil cells, based on its 
intense Raman TO band (1054 cm -1 at ambient con- 
ditions) which is well separated from the signal of the 
anvils. This has recently motivated investigations of the 
Raman spectra under extreme conditions of P and T— 

As pointed out by Holzapfel 6 -, the accurate determina- 
tion of the equation of state (EOS) of solids with very 
large bulk modulus is a stringent test for the ruby sen- 
sor calibration, a matter which has been under intensive 
debate recently^&^ i 10 ' 1 ^ 12 . For highly incompressible 
solids like c-BN, experiments in the megabar range are 
required in order to extract reliable values of the bulk 
modulus -Bo and its first pressure derivative B' . At am- 
bient temperature, the equation of state has been so far 
investigated up to 34 GPa in a He pressure medium 1 , 
to 66 GPa in N 2 13 , and 120 GPa in methanol-ethanol^. 
Although the maximum pressure in the latter study is 
relatively high, the small number of measurements, their 
limited accuracy, and the use of a pressure media prone to 
nonhydrostatic stress above ~ 20 GPa— reflects in poorly 
constrained and questionable values of Bq and B' . 

There is to date no available EOS data at simultaneous 
high pressure and temperature. By contrast, several the- 
oretical studie d 16 ' 17 ! 18 have reported thermal properties 
at ambient and high pressures. These works relied on the 
assumption that c-BN can be treated as a quasi-harmonic 



solid; this still needs to be validated by experiment. 

We present here careful and accurate measurements of 
the EOS up to 160 GPa at 295 K and up to 80 GPa 
at high temperatures (500-900 K). We used helium or 
neon (at high T) as pressure transmitting medium since 
they are known to provide the closest approximation to 
hydrostatic compression in the megabar range 9 ' 15 . The 
thermal expansion at ambient pressure has also been in- 
vestigated to 950 K. We show that our data can be rep- 
resented by a simple Mie-Griineisen model. This allows 
us to extract all the relevant thermodynamic parame- 
ters. The EOS is then used to determine the isothermal 
Griineisen parameter of the Raman TO band. Accord- 
ingly, the pressure-scale based on this Raman band is 
reformulated using physical constraints. 

This article is organized as follows: in section II, we 
present the experimental procedure; section III is devoted 
to the presentation of the room-temperature equation of 
state; the thermal expansion at room pressure follows in 
section IV; the presentation and analysis of P-V-T data 
at simultaneous HP-HT are given in section V; in section 
VI, the thermal dependence of the isothermal Griineisen 
parameter of the Raman TO band is investigated and 
a new formulation of the Raman pressure scale is then 
presented. Section VII finally gives concluding remarks. 



II. EXPERIMENTAL PROCEDURE 

The present experiments were conducted with single 
crystals of c-BN of size ranging from 3 to 15 /im. These 
were selected from a powder batch using their Raman sig- 
nal as a criterium for good crystallinity: crystals which 
presented intense Raman first-order TO and LO bands 
with lorentzian shapes and small bandwidths (~ 5 cm -1 ) 
were selected. A few crystals were then loaded into the 
experimental volume of a diamond anvil cell, along with 
ruby and SrB407t Sm 2+ as pressure sensors. Care was 
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taken in order to position the samples and the pressure 
sensors within a few micron distance to each other and at 
the center of the diamond culet. We used diamond anvils 
with fiat culets of 0.1 to 0.4 mm and rhenium gaskets. He- 
lium was chosen as pressure transmitting medium for the 
room temperature experiments above 50 GPa and neon 
was used otherwise. In the P — T range of the present 
experiments, these two pressure transmitting media are 
known to provide the best approximation to hydrostatic 
conditions*^. 

Membrane diamond anvil cells (MDAC) designed for 
high-temperature operation were used. The cells could 
be fitted as a whole inside a ring-shaped resistive heater. 
To achieve temperatures above 800 K, a smaller, addi- 
tional heater made of a resistive wire coiled around a 
ceramic tube, is positioned around the anvil-gasket as- 
sembly. The temperature of the heaters are regulated 
within 1 K using commercial devices. An isolated, K- 
type thermocouple is fixed by ceramic cement with its 
head in contact with the diamond anvil, close to the gas- 
ket. The ensemble is heated in air or in a Ar-H 2 (2%) 
reducing atmosphere. Numerous previous experiments 
have shown that the temperature measured by the ther- 
mocouple is within 5 K of the sample temperature^ ' 19 ' 20 . 

Pressure was determined using the pressure shift of the 
luminescence lines of ruby (at 295 K) or SrP^Oy: Sm 2+ . 
The pressures reported here are based on HolzapfePs 
2005 ruby scale 11 , hereafter denoted H2005. We 
also compared the results obtained with other ruby 
calibration s 6 ' 10 ' 11 ' 12 ' 21 and, as discussed below, H2005 
was found to provide excellent consistency between 
present and literature data. The calibration of the 
SrBztC^: Sm 2+ sensor— , initialy based on the ruby scale 
from Ref.l2ll. was also modified to match the H2005 scale. 
For the measurements above 100 GPa at 295 K, the ruby 
signal was too weak, thus the pressure was determined 
from the equation of state of 4 He (Ref. [23|, with proper 
correction for the H2005 ruby scale). The volume of 4 He 
was calculated using the reflections present on the same 
diffraction patterns as the c-BN sample and excellent 
agreement between this pressure determination and that 
from ruby was observed below 100 GPa. 

Angular-dispersive x-ray diffraction experiments were 
performed on bcamline ID27 of the European Syn- 
chrotron Radiation Facility (ESRF, Grenoble, France). 
The monochromatic beam (A = 0.3738 A) was focused 
to a w 7 x 10 /Ltm 2 spot. Diffracted x-rays were collected 
by a MAR345 image plate while the MDAC was continu- 
ously rotated about the 0-axis by ±20°. The images were 
integrated using the fit2D program 2 ^. Between 6 and 
9 single-crystal reflections could be observed depending 
on the sample, with a resolution up to « 0.83 A. All of 
them could be indexed in the zinc-blende structure (space 
group F43m) reported for this material 2 ^. The lattice 
parameter a was refined with the program UnitCell 26 , 
using the measured d-spacings of all observed reflections. 
The uncertainty on a was on average 5 x 10~ 4 A. 
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FIG. 1: (Color online) Volume per atom of c-BN as a func- 
tion of pressure at 295 K. The open and solid circles are 
present measurements on samples in a neon and helium pres- 
sure medium respectively. The estimated error bars are within 
the symbol sizes. Triangles: Ref. Q (in methanol/ethanol); 
dashed line: Ref. Q] (in helium); squares: Ref. [l3| (in nitro- 
gen). The solid line is the fit to the present data using the 
Vinet equation^Vb = 5.9062(6) A 3 /at, B = 395(2) GPa, 
Bq = 3.62(5)]. 



TABLE I: Comparison between values for the zero-pressure 
isothermal bulk modulus Bo and its first pressure derivative 
B' obtained by fitting the present data at ambient temper- 
ature to different EOS models: Vinel^, second-order Birch- 
Murnaghan (BM)^ and HolzapfePs AP2 form ?9 . The zero- 
pressure volume Fo(295 K) was fixed to 5.9062 A 3 /at. 



Model 


Vinet 


BM 


AP2 


Bo, GPa 


395(2) 


396(2) 


397(2) 


B' 


3.62(5) 


3.54(5) 


3.50(5) 



III. THE ROOM-TEMPERATURE EQUATION 
OF STATE 

The volume per atom V of c-BN at 295 K (V = a 3 /8) 
was measured up to 57.8 GPa in a neon pressure medium 
and up to 162 GPa in a helium pressure medium. Ex- 
cellent agreement is observed between the two data sets. 
The results are gathered in Fig. Q]and Table HVl As ex- 
pected, no phase transition was observed in this pressure 
range. 

To extract the values of the zero-pressure isothermal 
bulk modulus Bq and its first-pressure derivative B' 
from the present data, we considered three different phc- 
nomenological models of equation of state: the Vinet 
formal, the Birch-Murnaghan second-order form 2 ^ and 
HolzapfePs AP2 equation 2 ^. These relations were least- 
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TABLE II: Comparison between values of Bo (in GPa) and 
B' of c-BN and natural diamond (C) at 295 K obtained with 
different ruby scales from the literature. Bo and B are deter- 
mined by fitting the present data (c-BN) and Occelli et al.'s 
data (C) to the Vinet model^. In all fits, the values of the 
zero-pressure volume was fixed to 5.9062 A 3 /at (c-BN) and 
5.6733 A 3 /at (C). \ 2 indicates the goodness of fit. 







c-BN 




Diamond 




Ruby scale 


Bo 


B'o 


x 2 


Bo 


B'o 


x 2 


MXB1986 21 


397(3) 


2.75(6) 


14.2 


447(3) 


3.00(7) 


4.5 


AGSY1989i 


394(3) 


3.93(7) 


21 


440(3) 


4.38(8) 


6.4 


H2003' i 


390(2) 


3.49(5) 


12.1 


435(2) 


3.87(4) 


1.5 


DLM20042 


395(3) 


3.28(6) 


16.8 


440(3) 


3.66(7) 


4.6 


H2005ii 


395(2) 


3.62(2) 


12.8 


443(3) 


3.97(5) 


1.5 


CNSS2005^ 


387(3) 


3.64(7) 


18.2 


432(3) 


4.05(7) 


5.3 


DO2007— 


398(2) 


3.35(5) 


13.6 


444(2) 


3.72(5) 


2.6 



FIG. 2: (Color online) Comparison between experimental 
data at 295 K from this work (circles), Ref. Q (triangles), 
Ref. Qll (squares) and Ref. Q] (dashed line). The fit to the 
present data using the Vinet EOS§£ is used as the reference. 
The top graph shows the difference in pressure and the bot- 
tom one, the difference in volume multiplied by a factor 100. 

squares fitted to the data. The results are listed in Table|T] 
and the fit using the Vinet equation^ is plotted in Fig.Q] 
In the compression range probed here (0.77 < V/Vq < 1) 
the three models fit equally well the data and give very 
similar values for B (395(2) 397(2) GPa) and its first 
pressure-derivative B' Q (3.50(5) -^3.62(5)). Here as in all 
the text, the number in parenthesis indicate the standard 
deviation for the last digit of the fitting parameter; it 
does not reflect its absolute uncertainty, which primarily 
depends here on the one of the ruby standard. The vol- 
ume and pressure difference between experimental data 
and the Vinet fit is shown in Fig. [5) The volume Vq at 
P = and 295 K was taken as the measured value at 

1 atm of 5.9062(6) A 3 /at, corresponding to a lattice pa- 
rameter oo of 3.6152(2) A. The latter values are identical, 
within uncertainties, to previous reports 1 - ! 14 ' 30 . 

Figs. Q] and also shows experimental data obtained 
in previous studies 1 -^. For comparison with the present 
work, the original pressure values have been rescaled onto 
the H2005 scale. Aleksandrov ct al, 1 reported their mea- 
surements up to 34 GPa in He in the form of a quadratic 
fit of the compression factor Ap/ pn-where p is the den- 
sity, vs. the relative shift of the ruby R\ line. Their 
absolute volume values differ from ours by less than 

2 x 10~ 2 A 3 /at, corresponding to a maximum pressure 
difference of 2 GPa. By contrast, systematic deviations 
are observed between our measurements and those of 
Solozhenko et aL— and Knittle et aL— above about 40 
GPa, where their measured volume at a given pressure 
is systematically higher than our. This is likely due, at 
least in part, to the larger nonhydrostatic stresses pro- 



duced by the N2 or methanol-ethanol pressure medium 
used in Solozhenko et al^ and Knittle et al^ experi- 
ments respecitvely. 

As mentioned above, the values of Bo and B' Q obtained 
by fitting the EOS of highly incompressible solids like 
c-BN and diamond are very sensitive to the chosen pres- 
sure calibration for the ruby standard. As a matter of 
fact, the first hint that the commonly used Mao et al.'s 
1986 calibration^, (hereafter denoted as MXB1986) be- 
comes increasingly wrong with pressure was given by the 
measurements of the diamond EOS by Aleksandrov et 
alX to 40 GPa in 1989. More recently, the very accu- 
rate measurements of the diamond EOS up to megabar 
pressures by Occelli et aL— evidenced more firmly the dis- 
crepancy between the value found for B' using ultrasonic 
techniques (4.0 ± 0.7), on one hand, and that obtained 
from the EOS (3.0) when pressure is calculated with the 
MXB1986 scale, on the other hand. Consequently, a re- 
vision of the latter scale has been proposed by several 
authors^ ^ 10 ' 11 ! 12 , based on the diamond data and/or 
those recently obtained for the EOS of several metals 
that showed the same trend£. 

In Table [TTI we compare the Bo and B' parameters ob- 
tained by fitting the Vinet EOS to the present c-BN data 
using the various ruby calibrations proposed in the liter- 
ature. For comparison and consistency check, we also 
show the results obtained with the diamond data of Oc- 
celli et alZ. To decide which ruby scale(s) seem(s) more 
reasonable, we list in Table IIIII the values of Bq and B' 
determined by other means, whether experimental (ultra- 
sonic or Brillouin scattering measurements) or theoreti- 
cal. Unfortunately, there exists only one determination of 
B' for diamond using ultrasonic experiments under pres- 
sure, which has a limited accuracy, and none has been 
reported yet for c-BN. The largest source of compar- 
ison comes thus from the numerous theoretical studies 
based on first-principles (ab initio) techniques. Whereas 
the bulk modulus and zero-pressure volume show large 
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TABLE III: Literature data for the values of Vb, B and B' Q of c-BN and diamond ( 12 C). Vq is in A 3 /at and Bo in GPa. 
The various abbreviations are: DFT: density functional theory, LDA: Local density approximation, GGA: generalized gradient 
approximation, ZPE: zero-point energy, DG: Debye-Griineisen model, LD: lattice dynamics, VMC: variational quantum Monte- 
Carlo, DMC: diffusion quatum Monte-Carlo; "static" stands for calculations of the static lattice. 









c-BN 










12 C 




Reference 


V 


Bo 


B' 


Method 


Reference 


V 


Bo 


B'o 


Method 


Ref. 31 




400(20) 




Brillouin 


Ref. 32 




442(5) 


4.0(7) 


Ultrasound 


Ref. 16 


5.797 


395 


3.65 


DFT+LDA (static) 


Ref. 33 




444.8(8) 




Brillouin 




5.884 


387 


3.66 


DFT+LDA+ZPE (0 K) 


Ref. 34 


5.497 


473 


3.5 


DFT+LDA(static) 




5.888 


385 


3.66 


DFT+LDA+ZPE+DG(300 K) 


Ref. 8 


5.510 


465 


3.63(3) 


DFT+LDA (static) 


Ref. 18 


5.745 


391 




DFT+LDA+LD (300 K) 




5.697 


433(2) 


3.67(3) 


DFT+GGA (static) 


Ref. 35 


5.788 


397 


3.6 


DFT+LDA (static) 


Ref. 36 


5.529 


454 


3.65 


DFT+LDA+ZPE+QHA(300 K) 


Ref. 17 


5.742 


398 




DFT+LDA (static) 




5.722 


422 


3.72 


DFT+GGA+ZPE+QHA (300 K) 


Ref. 37 


5.718 


397 


3.59 


DFT+LDA (static) 




5.604 


472(4) 


3.8(1) 


VMC+ZPE+QHA (300 K) 


Ref. U 


5.954 


368 


3.6 


DFT+LDA (static) 




5.711 


437(3) 


3.7(1) 


DMC+ZPE+QHA (300 K) 


Ref. 38 


5.905 


370 


3.8 


Tight-binding 













variations among the various theoretical studies, the val- 
ues for B' appear much less dependent on the theoretical 
approaches and approximations used in each of them, es- 
pecially if we consider the most recent works. Indeed, 
the two latest studies of diamond give values of B' Q be- 
tween 3.65 and 3.8(1), whether the density functional 
theory (DFT), within LDA or GGA approximations, or 
the quantum Monte-Carlo approach is used. Similarly, 
all DFT-LDA calculations on c-BN give B' = 3.6 ± 0.05, 
irrespective of the considered pseudopotential. In the 
case of c-BN also, there is a nice agreement between the 
latest studies on the values of Bo which vary from 395 
to 398 GPa, i.e in very good agreement with the Bril- 
louin scattering experiment 3 ^ (400 ± 20 GPa). We see 
from Table |TT] that the ruby scales H2003, H2005 and 
DO2007 are those giving the best consistency with avail- 
able literature data, while providing the best fits based 
on the least-squares \ 2 criterium. The DLM2004 and 
CNSS2005 calibrations give close results but a larger % 2 , 
which could be ascribed to the used functional form. 
This also confirms that the MXB1986 scale becomes in- 
creasingly wrong with pressure, underestimating it by 
about 8.8% at 160 GPa. We note that H2003, H2005, 
DLM2004, CNSS2005 and DO2007 scales agree within 
3% at 160 GPa, which may be considered as the present 
uncertainty of the ruby scale at this pressure. It is 
now clear that independent experimental determination 
of B' , such as given by sound-propagating experiments, 
for both c-BN and diamond could help to better establish 
the ruby scale in the megabar range. 



IV. THERMAL EXPANSION AT AMBIENT 
PRESSURE 

We performed volume measurements at room pressure 
as a function of temperature between 295 and 948 K. The 
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FIG. 3: (Color online) Volume and thermal expansion co- 
efficient as a function of temperature. The solid circles are 
present data and triangles show measurements of Ref.l3(il. The 
fit to both data sets using Eq. ([1} is shown as the red line. 
The blue dotted line represents the calculations of Albeit. In 
the inset, the volume thermal coefficient ceo is represented. 



results are plotted in Fig. [3] and reproduced in Table HVl 
They are compared to those of Slack and Bertram 3 ^ be- 
tween 77 and 1289 K. The two data sets agree within 
error bars in the overlapping range, although we system- 
atically find larger volume values at a given temperature 
above 300 K. The calculations of Albe^ are also shown as 
the dotted line in Fig. (3J They are based on the density 
functional theory (DFT) within the local density approxi- 
mation (LDA) for the static part, and a Debye-Griineisen 
model for the thermal part. The calculated volume fol- 
lows very well the experimental data below 500 K but 
increasingly overestimate them above this temperature. 
Present and Slack and Bertram?^ data were fitted to- 
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TABLE IV: Experimental P-V -T data of c-BN obtained 
in the present work. P is in GPa, a in A, V in A 3 /at and T 
in K. 



p 


T 


a 


V 


P 


T 


a 


V 


0.0 


295 


3.6152 


5.907 


0.0 


377 


3.6158 


5.909 


1.5 


295 


3.6110 


5.886 


0.0 


479 


3.6174 


5.917 


4.5 


295 


3.6029 


5.846 


9.8 


497 


3.5879 


5.773 


7.4 


295 


3.5938 


5.802 


21.6 


499 


3.5555 


5.618 


9.3 


295 


3.5880 


5.774 


2.0 


500 


3.6110 


5.886 


11.1 


295 


3.5832 


5.751 


25.1 


500 


3.5465 


5.576 


12.7 


295 


3.5792 


5.732 


13.9 


500 


3.5760 


5.716 


13.0 


295 


3.5793 


5.732 


46.5 


501 


3.4983 


5.352 


13.5 


295 


3.5771 


5.721 


44.5 


501 


3.5024 


5.370 


15.4 


295 


3.5715 


5.695 


0.0 


573 


3.6181 


5.920 


17.8 


295 


3.5658 


5.667 


39.3 


600 


3.5145 


5.426 


19.9 


295 


3.5607 


5.643 


47.9 


600 


3.4965 


5.343 


21.4 


295 


3.5578 


5.629 


56.5 


600 


3.4791 


5.264 


22.3 


295 


3.5541 


5.612 


66.3 


600 


3.4609 


5.182 


24.4 


295 


3.5494 


5.590 


76.2 


600 


3.4458 


5.114 


25.0 


295 


3.5484 


5.585 


84.2 


600 


3.4264 


5.028 


26.0 


295 


3.5461 


5.574 


0.0 


674 


3.6201 


5.930 


28.6 


295 


3.5395 


5.543 


8.2 


748 


3.5962 


5.814 


38.6 


295 


3.5164 


5.435 


4.6 


748 


3.6069 


5.866 


42.7 


295 


3.5040 


5.378 


18.9 


749 


3.5653 


5.665 


48.2 


295 


3.4930 


5.327 


35.8 


749 


3.5243 


5.472 


51.9 


295 


3.4849 


5.290 


42.1 


750 


3.5095 


5.403 


53.1 


295 


3.4823 


5.278 


49.5 


750 


3.4938 


5.331 


57.8 


295 


3.4746 


5.244 


14.7 


750 


3.5789 


5.730 


66.6 


295 


3.4565 


5.162 


0.0 


775 


3.6216 


5.938 


77.3 


295 


3.4405 


5.091 


0.0 


873 


3.6233 


5.946 


87.0 


295 


3.4249 


5.022 


54.2 


900 


3.4841 


5.287 


95.7 


295 


3.4087 


4.951 


0.0 


948 


3.6246 


5.952 


101.6 


295 


3.3989 


4.908 










109.8 


295 


3.3860 


4.853 










119.3 


295 


3.3722 


4.793 










119.7 


295 


3.3721 


4.793 










124.3 


295 


3.3634 


4.756 










133.3 


295 


3.3525 


4.710 










139.7 


295 


3.3437 


4.673 










147.7 


295 


3.3334 


4.630 










154.5 


295 


3.3246 


4.593 










162.5 


295 


3.3147 


4.553 











gether using the second-order approximation to the zero- 
pressure Griineisen equation of state (see Ref. [H and 
references therein). In this approximation, the tempera- 
ture dependence of the volume is given by: 



Vo(T) - Vb.o 



U th (T) 



(1) 



Q-bU th (T) 

where Uth(T) is the internal energy due to lattice vi- 



brations, b = \(B' QQ - 1) and Q = (V , B ofi )/j. Vb,o, 
B o and B' are respectively the volume, bulk modulus 
and its first derivative at zero pressure and temperature. 
7 is the thermobaric Griineisen parameter defined by: 



7 = -(V/U th )(dF th /dV) T = (V/U th )P th (2) 

where Fth(V,T) is the thermal part of the Helmholtz 
free energy and P t h = — (dF t h/dV)T- If 7 is inde- 
pendent of T, which we assume here, then 7 = 7^ 
where 7^ is the thermodynamic Griineisen parameter 
(ith = (V/C v )(dP/dT)v with C v the specific heat at 
constant volume) 39 . Uth{T) was evaluated within the 
Debye model: 



e /T 



exp(x) — 1 



dx (3) 



where R is the ideal gas constant and Oo the Debye 
temperature at P = 0. For the latter, we used the 
value of 1700 K estimated from the infrared spectrum 
by Gielisse et al. 40 , which is also in agreement with data 
on specific hea t 18 ' 41 . i?o,o was fixed at 397 GPa, which 
is obtained by adding to the value of -B (295 K) found 
above (Vinet fit), the small correction (2 GPa) calcu- 
lated by Albeit between and 300 K (Table [III]). B' 
was fixed equal to B' (295 K) = 3.62. 

The fit of the V (T) data with Eq. {TJ) then gives Vo,o — 
5.9026(4) A 3 /at and ltm =1.04(1) ( 7tM = lth(P = 0)). 
It is shown as the solid red line in Fig. [3l The tempera- 
ture dependence of the thermal expansion coefficient ao, 
calculated by differentiating Eq. |T]), is plotted in the in- 
set of the same figure. It is seen that Eq. (fTJ) gives a 
good reproduction of Vq(T) in the considered tempera- 
ture range. We also note that the present value for j th0 is 
close to the one obtained by Brillouin-zone integration of 
the phonon dispersion curve (< 7„ >= 0.95) determined 
by Kern et al¥- using density functional theory. 



V. HIGH PRESSURE AND TEMPERATURE 
EQUATION OF STATE 

High pressure and temperature measurements were 
performed up to 80 GPa at 500, 600, 750 and 900 K, 
with the results listed in Table HVl These are to our 
knowledge the first reported data at simultaneous HP- 
HT conditions. The P-V-T data set is represented in 
Fig. Ha). 

This data set was least-squares fitted using a Mie- 
Griineisen model. In this model, the pressure is expressed 
as the sum of a "cold" (Pq) and thermal (Pth) part, i.e.: 



P(V,T) = P (V) + P th (V,T) 



(4) 



with P th (V, T) = at K. To represent P (V), a Vinet 
equation 27 is used, whereas the thermal part is written 
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FIG. 4: (Color online) (a) Experimental volume of c-BN as a 
function of pressure and temperature obtained in the present 
work. Crosses show experimental data and the dashed-line 
is the Vinet fit to the data at 295 K. (b) Difference be- 
tween experimental pressures and predicted ones using the 
Mie-Griineisen model with the parameters given in Table FVl 



in the Debye-Griineisen (quasi-harmonic) approximation, 
using Eqs. (J2]t3j) . In this approximation, we have the 
identity jth = 7 = 7d, where jd = — <9 ln.(0)/9 ln(V) 
is the Debye-Griineisen parameter. The variation of 
the Debye temperature with volume is thus given by 
that of "fth- Here jth was allowed to vary with the 
compression ratio according to the empirical relation 
7th = Jtho(V/Vo,o) 9 , where q was taken as constant^. 
In the fit, the only varying parameter is q. All the others 
(V^o, B q, B' q , 7 t ^o and O ) were fixed to the values 
determined as described above. The obtained value for q 
is 4 ±1.5 and the whole parameter set is given in Table fVl 
Fig. 0Jb) shows the difference between the experimental 
pressures and the one predicted by the present model. 
This difference is for most data points smaller than ±1 
GPa and the rms deviation is 0.6 GPa. If we allow i?o,o 
and jtho to vary in the fitting procedure, the resulting 
values are identical to the starting ones within their stan- 
dard deviations. The parameter set given in Table fVl can 
thus be regarded as the best one for the Mie-Griineisen 



TABLE V: Parameters for the P-V-T equation of state 
of c-BN based on a Mie-Griineisen model, with a Vinet 
equation^ for the static part and a Debye-Griineisen model 
for the thermal part [Eq. Q]. 



T/ ,o,A 3 /at 


Bo,o, GPa Bq 


e 0> k 


itho q 


5.9026(4) 


397(2) 3.62(5) 


1700 


1.04(2) 4(1.5) 
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FIG. 5: Variation of the thermal expansion coefficient with 
temperature along several isobars, as deduced from the Mie- 
Griineisen EOS. 



EOS based on the present data. 

The value found for q is large compared to the typical 
range for this parameter (0.8 < q < 2.2)^2. This is to be 
related to the rapid decrease of the thermal expansion 
coefficient with pressure that is observed from present 
data. Indeed, jth an d a are related by the thermody- 
namic relation a = (l/B)"f t h,Cv- In the quasi-harmonic 
approximation, 7^ may be obtained by the Brillouin- 
zone integration of the phonon mode Griineisen param- 
eters. A large reduction of a may thus indicate that the 
Griineisen parameters for the acoustic modes, which give 
the dominant contribution, rapidly decrease with pres- 
sure. As a matter of fact, the calculated 7 for the trans- 
verse acoustic branches at P — are very low (~ 0.25) 
along some directions of the Brillouin zone, especially 
near the zone boundaries^. In diamond, which presents 
similar phonon dispersion curves, calculations show a de- 
crease of these parameters with pressure, eventually lead- 
ing to negative values of 7 and a at ultra-high pressures 
(P > 700 GPa) 43 . We also note that a similar phe- 
nomenum is responsible for the negative thermal expan- 
sion in Si at ambient pressure and low temperatures^. 
To illustrate the decrease of a with pressure, we plot in 
Fig. [5] its variation with temperature along several iso- 
bars from to 100 GPa as deduced from the present 
Mie-Griineisen EOS. 
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VI. THE CONSTRAINED RAMAN PRESSURE 
SCALE 

A. Temperature dependence of the TO mode 
Griineisen parameter 

The present construction of a P-V-T EOS enables us 
to reanalyze our previous measurements^ of the Raman 
TO mode at HP-HT. In particular, we are now able to 
precisely determine the isothermal mode Griineisen pa- 



rameter 7 io (T) 



where v TO is 

the frequency of the TO mode. To do so, the sample 
volume was calculated for each P-T conditions at which 
v TO was measured (see Ref. H) by inverting Eq. @ . 

The results are reported in Fig. [SI where we plot the ex- 
perimental values of ln(V TO ) as a function of ln[V/Vb(T)]. 
In the scanned pressure range (P < 21 GPa), and for each 
of the five studied isotherms (300 < T < 723 K), these 
two quantities appear to be linearly related. This tells 
us that 7 (T) is constant along each isotherm and may 
then be directly determined by a linear regression of the 
data, using the expression: 



ln*/ TO (T/T) = -7 TU (T)ln 



V(P, T) 

w 



+ lnz, TO (T) (5) 



where v$°{T) 



TO 



(P = 0, T). The fits to the latter 



equation are shown as dotted lines in Fig.[S] The values of 
7 TO (T) and vfi°(T) so obtained are listed in TableEE] It 
can be seen that, within error bars, 7 TO (T) is constant in 
this P-T range, with an average value of 1.257(5). This 
value is slightly larger than that given by Aleksandrov et 
alX [1.188] and in good agreement with the theoretical 
determination of Kern et al*Z- [1.2]. 

It can also be observed that Vq°(T) follows Vq(T) in 
a linear way in the probed T range. Extrapolating this 
line to K gives a value of Vq° at zero temperature of 
1055.1(1) or 1054.6(2) cm -1 , depending on whether the 
fitted or experimental Vq°(T) is used. 



B. The high-temperature pressure scale 
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FIG. 6: (Color online) Frequencies of the TO Raman mode 
as a function of V/Vo(T) on a logarithmic scale. The crosses 
show experimental data£ and the dashed lines are linear fits. 



TABLE VI: Isothermal mode Griineisen parameter of the Ra- 
man TO mode based on the Raman data reported in Ref. 
and present P-V-T data. Vo is the volume at P = obtained 



by inverting Eq. (j4| 



The value of Vq° 



to 



(P = 0, T) pre- 



dicted by the linear regression of the data [Eq. [5] is also given 
and compared to the direct measurements. T is in K, Vo in 
A 3 /at and Uq° in cn 
by linear regression of u£°(T) vs. V (T). 



The values at T — are obtained 



T 


Vo 


v™ (fit) 


vl° (exp.) 


1 T ° (fit) 


300 


5.9055 


1054.0(1) 


1053.93(10) 


1.251(3) 


373 


5.9084 


1052.9(2) 


1052.43(20) 


1.256(5) 


473 


5.9139 


1050.8(2) 


1050.17(13) 


1.262(6) 


573 


5.9209 


1048.3(1) 


1048.17(27) 


1.259(2) 


673 


5.9290 


1045.1(1) 


1045.45(13) 


1.262(4) 


723 


5.9334 


1043.6(1) 


1043.91(41) 


1.253(4) 





5.9026 


1055.1(1) 


1054.6(2) 





In Ref. [U we recalled the various reasons that make 
c-BN a good candidate for pressure measurement in a 
DAC at high temperature. A pressure scale was given 
based on the measurements of the TO mode frequency 
(v TO ) at high P-T . The form of this pressure scale de- 
rives from the first-order Murnaghan EOS, and is ob- 
tained by inverting the following equation for the P-T 
dependance of v TO : 

^o (P)T) = ^o (T) ^ 1 + _^_ p j (6) 

In Ref. S, B (T) and B' were considered as fit param- 
eters, taking for 7™ the value reported by Ref. U We 



assumed then that 7 was independent of temperature, 
which is confirmed by the present work. Our EOS allows 
us to better constrain Eq. ([6]) by imposing physical con- 
straints on the various parameters. Since this pressure 
scale is to be used for high temperatures, we only con- 
sider variations above 300 K. The variation of Bq with 
temperature was deduced from our Mie- Griineisen EOS 
by fitting isotherms at 100 K intervals by a Vinet equa- 
tion. A quadratic form is found suitable to represent the 
behavior of Bo between 300 and 2000 K, with the follow- 
ing expression: 

S (T> 300K) = 396.5(5)- 0.0288(14) (T- 300) (7) 
-6.84(77)10- 6 (T-300) 2 
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FIG. 7: (Color online) Comparison between the present pres- 
sure scale [Eq. ©] and the one given by Goncharov et al. 
[Ref.H Eq. (2)]. 



The data for v TO {P,T) was then fitted to Eq. ©, us- 
ing Eq. as the expression for B (T). Adding a tem- 
perature dependence to B' does not improve the fit, so we 
kept it as constant. The value of 7™ = 1.257 obtained 
above was used. AsinRef.i, v£°{T) was expressed as a 
quadratic form, where the coefficients are allowed to vary 
within the boundaries given by Ref. |45| . We thus obtain 
a revised pressure scale, reading as: 



P = (fl (T)/3.62) 



TO 



(P,T) 



2.S70 



- 1 



(8) 



with B (T) as in Eq. Q and i#°(T) = 1058.4(5) - 
0.0091(23) T - 1.54(22) x 10~ 5 T 2 

Recently, Goncharov al. 5 reported Raman measure- 
ments up to 1750 K and 40 GPa in argon pressure 
medium and proposed a pressure scale in a form simi- 
lar to Eq. IS} . The two scales are compared in Fig. [7| in 
the pressure range 0-50 GPa at 300, 1000 and 1500 K. 
The difference increases with pressure and temperature 
and reaches ~10% at 50 GPa and 1000 K. The fact that 
the present scale is constrained by measured physical pa- 
rameters gives us confidence that it can be safely used 
in the 0-100 GPa and 300-1000 K range. P-V-T as well 
as Raman data to higher T would be desirable to extend 
the calibration. 



VII. CONCLUSIONS 

We have reported experimental P-V-T data on c-BN 
from x-ray diffraction experiments in a resistively-heated 
diamond anvil cell. Volume measurements extend to 160 
GPa at 295 K and 80 GPa at 500-900 K. To our knowl- 
edge, these are the first reported EOS data at simultane- 
ous HP-HT conditions. By fitting the room-temperature 
EOS to various EOS models, we extracted the values 
of the bulk modulus [395(2) GPa] and its first-pressure 
derivative [3.62(5)]. We have also examined the varia- 
tion of these parameters with respect to the chosen cal- 
ibration of the ruby standard, and showed that the one 
recently proposed by Holzapfelii (H2005) provides very 
good consistency between static compression data and 
ab initio predictions. Coupling the present information 
with independent measurements of B' , such as obtained 
by sound-propagating experiments, would be of great in- 
terest to better constrain the ruby standard calibration in 
the megabar range. A good description of our full data 
set was obtained using a Mie-Griineisen model, where 
the thermal pressure originating from lattice vibrational 
energy is calculated in the Debye approximation. We ob- 
serve a rapid decrease of the thermal expansion coefficient 
with pressure, which itself reflects on the strong variation 
of the thermal Gruneisen parameter. The Mie-Griineisen 
EOS was then used to determine the mode Gruneisen pa- 
rameter of the Raman TO mode, which was found to be 
temperature independent in the range 300-723 K. A new 
formulation of the Raman pressure scale was then de- 
duced from the present results, which should hold in the 
P-T range 0-100 GPa and 300-1000 K. Extension of this 
work to higher temperatures would be valuable to better 
constrain the thermal effects, and extend the calibration 
of the Raman scale. 
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